Thermoelectric performance of granular semiconductors 
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We study thermoelectric properties of granular semiconductors with weak tunneling conductance 
between the grains, gt < 1. We calculate the thermopower and figure of merit taking into account 
the shift of the chemical potential and the asymmetry of the density of states in the vicinity of the 
Fermi surface due to n- or p-type doping in the Efros-Shklovskii regime for temperatures less than 
the charging energy. We show that for weakly coupled semiconducting grains the figure of merit 
is optimized for grain sizes of order 5nm for typical materials and its values can be larger than 
one. We also study the case of compensated granular semiconductors and show that in this case 
the thermopower can be still finite, although two to three orders of magnitude smaller than in the 
uncompensated regime. 
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A major research area in nano-science and materials 
physics is the search for highly efficient thermoelectric 
materials. Although there has been extensive research ef- 
forts in the last several decades, the progress in this quest 
has been limited until recently. It was found that for fur- 
ther improvement in efficiency inhomogeneous/ 'granular 
thermoelectric semiconductors are especially suited 
These materials are now accessible as next generation 
devices for conversion of thermal to electrical energy and 
vice versa and technologically important due to the pos- 
sibility of direct control of the system parameters. The 
dimensionlcss figure of merit, ZT = S 2 <jT/k, is the 
preferred measure for the performance or efficiency of 
thermoelectric materials, where S is the thermopower 
and a and k the electric and thermal conductivities, re- 
spectively 0, El- Recently, ZT values of 2.4 in layered 
nanoscale structures 0] at 300K, and 3.2 for a bulk semi- 
conductors with nanoscale inclusions [B| at about 
were reported. These high values of ZT are in the range 
for applications and therefore call for the development 
of a comprehensive quantitative description of thermo- 
electric properties of granular semiconductors, which can 
serve as basis for a new generation of thermoelectrics. 

In this paper we investigate the thermopower S and the 
figure of merit ZT of granular semiconductors focusing on 
the case of weak coupling between the grains, gt < 1, see 
Fig. [T] Each semiconducting nanocrystal is characterized 
by two energy scales: (i) the mean energy level spacing 
S = l/(na d ), where n is the density of states at the Fermi 
surface, a is the grain size, and d is the dimensionality of 
a grain, and (ii) the charging energy E c = e 2 /(e r a) with 
e r being the dielectric constant. In semiconductors the 
density of states n is of about two orders of magnitude 
smaller than that in metals. Thus in semiconducting 
dots 5 can be of order of the charging energy, S ~ E c , 
in contrast to metallic granular materials where typically 
8 <C E c 6] . Our considerations are valid for temperatures 
T <E C . 

The internal conductance of a grain is taken larger than 
the inter-grain tunneling conductance, gt, which controls 



macroscopic transport properties of the sample Q. In 
this paper we consider g t < 1, i.e., smaller than the 
quantum conductance, which is the typical experimen- 
tal situation 0, 0] . 




FIG. 1: (Color online) Sketch of a nanogranular material 
showing typical electron (e) and hole (h) transport. In the 
upper row of grains an inelastic electron tunneling process is 
shown and in the lower row a co-tunneling loop responsible 
for the electronic part of the heat transport is presented. The 
energy (e) transport goes from the "hot" (H) to the "cold" 
(C) side. 

In the case of diagonal (short range) Coulomb in- 
teraction, the total probability for an electron tunnel- 
ing through many grains via elastic or inelastic co- 
tunneling can be written as the product, P = fl^Li -Rj 
of the individual probabilities of single elastic/inelastic 
co-tunneling events through single grains with N = r/a 
is the number of grains. The probability P is related to 
the localization length £ as P ~ e~ r l^ . 

Semiconducting nanocrystal arrays are described by 
the Hamiltonian 
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where i,j are the grain indexes and the summation in 
the second term of the r.h.s. of Eq. |T]) is performed over 
nearest neighbors. The term Tiy< is the Hamiltonian for 
the single grain i including the free electron energy and 
the diagonal Coulomb interaction, and the second term 
is the tunneling Hamiltonian between the adjacent grains 
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i and j with Uj being random tunneling matrix elements 

and cl [cj.*?] the creation [annihilation] operator on the 
ith grain. Due to the large mean energy level spacing in 
semiconducting grains S ~ E c , only a few terms of the 
fc-sums are important. 

In Ref. [f| it was shown that the probability for elastic 
Pf and inelastic Pj n co-tunneling through an array of 
weakly coupled semiconducting grains has the form 



p? 



1 



[rmax((5,S c )] 2 ' 
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(2) 



where r is the electron escape time from a grain. Thus, 
for the elastic/inelastic localization length we obtain 



i 



a/2 



ln[rmax(<5, E c 



a/2 



ln[Tmax(6,E c )] +S/T ' 
(3) 

The localization length £ cl / m j s related to the characteris- 
tic temperature scale To = e 2 / (e r ^ el ' ln ), which is of order 
EclnP" 1 . Below we derive the thermopowcr, S, ther- 
mal conductivity, n, and figure of merit, ZT for granular 
semiconductors using Eq. ([3]). 

One remark is in order: The thermopowcr S of lightly 
doped compensated semiconductors was investigated in 
the past 0,0, 11 , 12 . However, all previous studies were 
concentrated on the Mott variable range hopping (VRH) 
regime, with conductivity being <r(T) ~ exp (Tm/T) 1 / 4 , 
where Tm is the Mott temperature [l3j . In granular ma- 
terials the Mott VRH regime is hard to observe. Indeed, 
in semiconductors the Efros-Shklovskii (ES) law 14, 15| 




FIG. 2: (Color online) Simulation results for the density of 
state n(e) vs. energy for different filling factors v, indicated 
by the arrow across the curves from i/ = 0.5to^ = 0.8in steps 
of Av = 0.05 (u > 1/2 corresponds to n-type doped semicon- 
ductors). The simulations were done for a 2D Coulomb glass 
model, simulating the whole weakly coupled granular sam- 
ple Hi, E3 for a system size of 500 2 . The inset shows the 
dependence of the shift of the chemical potential A/i and the 
asymmetry of the density of states An on v (see text). Using 
these data we extract the numerical coefficients a\ m 2.9 and 
ct2 ~ 0.4 in Eq. by linear fitting. 



may turn into the Mott behavior with an increase of tem- 
perature. This happens when the typical electron energy 
e involved in a hopping process becomes larger than the 
width of the Coulomb gap A c , i.e. when it falls into the 
flat region of the density of states where Mott behavior 
is expected. To estimate the width of the Coulomb gap 
A c , one compares the ES expression for the density of 
states n(A c ) ~ (e r /e 2 ) d |A c | d_1 with the DOS in the ab- 
sence of the long-range part of the Coulomb interaction, 

no (d — 2, 3 is the dimensionality of a sample). Using the 

i 

condition n(A c ) ~ no we obtain A c = (noe 2d /ef) d_1 . In- 
serting the value for the bare DOS, no = l/E^^, into the 
last expression we finally obtain A c ~ E c . This means 
that there is no flat region in the DOS for T < E c . 

Here we discuss two effects: i) we calculate the ther- 
mopower S of granular semiconductors taking into ac- 
count the shift of the chemical potential A// = a\(y — 
l/2)To with v being the electron filling factor {y is re- 
lated to the compensation level of semiconductors, here 
we concentrate on n-type doped semiconductors with 
1/2 < \v\ < 1), ai a dimensionless numerical coeffi- 
cient, and the asymmetry of the density of states (DOS) 
An = «2 [y— 1/2)T " 1 with a-i being a numerical constant; 
ii) we show that even in the absence of the chemical po- 
tential shift, Au = 0, and asymmetry of DOS, An = 0, 
the thermopowcr S is still finite, although small, due to 
co-tunneling processes. We start with the former case. 

To calculate the thermopower of granular materials in 
the regime of weak coupling between the grains it is nec- 
essary to take into account electrons and holes because 
the contributions of electrons and holes cancel in the lead- 
ing order. In general the thermopower is proportional to 
the average ener gy t ransferred by charge carriers and can 
be written as [{jl lift 11, 12 1 



S 



1 

2eT 



[(e-A)e + {e-p)h] 



(4) 



Here the subscripts e and h refer to electrons and holes 
and jl = fi + Au is the shifted chemical potential. The 
expression in the square brackets of the r.h.s. of Eq. ((4]) 
describes the average energy transferred by charge carri- 
ers (electron or hole) measured with respect to the shifted 
chemical potential jl. The average energy in Eq. (j4j can 
be calculated as follows 



{e - fj) e /h = — 



J de (e - pj)n(e - p)f B / h {e - p)e 



J den(e - jl)f e / h {e - p)t 



2A^ 



(5) 

Here / e /ft,( £ ) is the Fermi function for electrons or holes, 
A = \/TqT the typical transfer energy in one hop, and 
n{e) is energy dependent the DOS. As we will see later, 
it is crucial to take into account the asymmetry of the 
DOS and the shift of the chemical potential in order to 
obtain a finite result in Eq. (fj]) . 
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The DOS, n(s) in Eq. for p, — A c < e < /i has the 
following form [X 01 ] 



n(e) oc |e- Mr" 1 ! 1 ~ (e ~ m) a ™] 



(6) 



and is constant, no, outside the Coulomb gap region, 
e < p, — A c and e > ft, where A c ~ E c is the width 
of the Coulomb gap. The shift of the chemical potential 
A/i = p, — fj, and the asymmetry of the DOS, An, are 
explicitly defined above Eq. (j4|). 

To support our choice for the expression of the DOS 
n(e) in Eq. ([6|) we numerically compute the DOS for a 
2D Coulomb glass model to simulate the whole system of 
semiconducting grains (see Refs. 0, 17 1 for details) at 



arbitrary filling factor v using first principles. The result 
of the simulations is shown in Fig. [2] These simulations 
clearly indicate that for a filling factor v ^ 1/2, the DOS 
is asymmetric and the chemical potential is shifted. Us- 
ing these results we can identify the dimensionless nu- 
merical coefficients a\ ~ 2.9 and ai ~ 0.4 by a simple 
linear fit. We note that ai > a2, thus the contribution 
to the thermopowcr caused by A/i in Eq. is dominant. 

Now, we can calculate Eq. ([5]) and the analog contri- 
bution for holes using Eq. ((6]) for the DOS together with 
the numerically found values for a\.i [20I ]. Finally wc 
derive the expression for the thermopower S of granular 
semiconductors in the limit of weak coupling between the 
grains 



d [ Am + AnToT] _ 1/2 -v 
eT e 



d 



T 
T 



a 2 



(7) 



We note that the r.h.s. of Eq. ((7|) vanishes for filling 
factor v = 1/2, i.e., for compensated semiconductors. 

Equation (J7]) is valid for temperatures T < T = 
-EclnP" 1 and weak coupling between the grains g t < 1. 
Under this condition, the electric conductivity a is @, 13 



2e 2 a 2 ' d g t e-V^. 



(8) 



The thermal conductivity k consists of two parts: the 
electron, n e and phonon, k p %. The phonon contribution 
K p h at temperatures T < Op, where Qjj is the Dcbyc 
temperature is given by jlil . flpj ] 



i 2 P h d T [t/q d 



d-1 



(9) 



where l p h = \f cxp(©£>/[aT/AF]) is the phonon mean 
free path in granular semiconductors with A^ being the 
Fermi length. [For a ~ lOnm, X F ~ lA, ® D — A50K, 
one obtains l p h — lnm at T ~ 200A'] The main contri- 
bution to the electric part n e of the thermal conductivity 
k appears due to a single closed co-tunncling loop (see 
Fig.[T]). An electron executing a co-tunncling loop brings 
back its charge to the starting grain and hence there is 
little change in the electrical conductivity and therefore 
the classical activation term is absent. However, there 




FIG. 3: (Color online) Diagrams describing lowest (a) and 
higher (b) order co-tunneling processes (multiple tunneling 
events are indicated by the dotted lines) . The solid lines de- 
note the propagator of electrons. The tunneling vertices are 
described by the circles. The wavy lines indicate the external 
coupling to the heat vertices. 



is no requirement that the returning electron has exactly 
the same energy (due to inelastic processes). The leading 
contribution to K e is proportional to g 2 and is depicted 
in the diagram shown in Fig. [5Ji). The analytical result 
corresponding to this process can be estimated as follows 



2 2- 
9t a 



'T 



-2S/T 



t max(<5, Ac)] 2 



(10) 



where we used Eq . (J2J for the inelastic co-tunncling prob- 
ability P™. Since at relatively high temperatures the 
phonon contribution k p h to thermal conductivity is dom- 
inant, we can neglect the contribution n e in the following. 

Substituting Eqs. fTJ) and ([9]) into the expression for 
the figure of merit ZT ~ S 2 aT '/ K p h we obtain the result 



ZT 



2d 2 g t e-V^ [A/jl/T + AnT ] 
(A i r/a) 2 - d (T/6 D ) d - 1 



,(d-2)- 



(11) 

Here the expressions for Au and An are given above 
Eq. Using Eq. (fTTj) we can calculate the tem- 

perature T* at which ZT has its maximum value, 
given by the solution of the quadratic equation T* = 
± (a+ (d+ 1)T*) 2 , where a = {2-d)Q D \ F /a. In d = 2 
we get — To/[4(d+ 1)], while in 3D the existence of 
a maximum depends on the values of a and To- 

In Fig. [3] we plotted the figure of merit ZT for a two- 
and three-dimensional system, using typical parameters 
for granular semiconductors: \ F ~ 1A, Op ~ 450K, 
T ~ 2E C , e r ~ 4, and g t ~ 0.1. Figure QJi) shows the 
temperature dependence of ZT for two different filling 
factors v = 0.6, 0.65 for a grain size of a ~ 5nm (For this 
size and the above dielectric constant we get E c ~ 800K). 
We clearly see that the figure of merit can well exceed 
one. However, we remark that ZT depends inversely 
proportional on the numerical coefficient of K p h, which 
is assumed to be one here. Figure [4)d) shows the grain 
size dependence of ZT at fixed temperature T = 200K. 
At this temperature for the above system parameters, 
the figure of merit is optimal for grain sizes a^ D ~ 4nm 



in 3D and a. 



2D 



6nm in 2D. 
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2D, v=0.65 




FIG. 4: (Color online) Plots of the figure of merit ZT vs 
temperature T (a) and grain size a (b) for two (2D) and 
three (3D) dimensional samples and two different filling fac- 
tors v — 0.6,0.65 (n-type doped granular semiconductors). 
Typical material parameters are given in the text. In a), the 
grain size is chosen to be a = 5nm and in b), the temper- 
ature is fixed at 200K. One sees that ZT can be optimized 
by choosing appropriate grain sizes: For the used parameters, 
these sizes are about 4nm in 3D and 6nm in 2D. 



Now, we concentrate on the compensated regime (fill- 
ing factor v = 1/2). In this case Eq. ([7]) predicts zero 
thcrmopower S. Thus, the fundamental question exists 
what mechanism may lead to a finite thermopower in this 
case. We start our consideration with the fact that the 
thcrmopower S can be expressed in terms of the ther- 
moelectric coefficient rj and the electric conductivity a as 
S = rj/a. Thus, to calculate the thermopower S one has 
to know the thermoelectric coefficient r\. To estimate f) 
we use the diagram shown in Fig. [2t>) (here the tilde in- 
dicates the compensated case). Since the dominant con- 
tribution to fj vanishes due to particle-hole symmetry, to 
obtain a nonzero result it is necessary to take into ac- 
count the fact that the tunneling matrix elements t and 
the density of states n depend on energy. In the leading 
order, the corrections to both quantities are proportional 
to T/Ep, where Ep is Fermi energy Q. As a result, the 
thcrmoelccric coefficient is given by the expression 



V 



,2-d 



9t [T/E F ] 



-y/To/T 



(12) 



where the temperature scale To was defined below 
Eq. ©. Substituting Eq. © and flU into the expression 
for thermopower we obtain 

S~(l/e){T/E F ). (13) 

It follows that the thcrmopower is finite although small, 



leading to a small figure of merit as well, since ZT ~ S 2 . 
Using Eqs. (J7J) and l|13p one can see that the ratio of 
two thermopowers for compensated [v = 1/2, Eq. I|13[) ] 
and for n-typc (v > 1/2, Eq. (|7|) regimes is of order 
S/S ~ T 2 /(T Q E F ) ~ 1(T 3 - 1CT 2 <c 1. 

Iii conclusion, we studied thermoelectric properties of 
granular semiconductors at weak tunneling conductance 
between the grains, gt < 1. We calculated the ther- 
mopower S and figure of merit ZT taking into account 
the shift of the chemical potential and asymmetry of the 
density of states for the exemplary case of n-typc doping. 
We showed that the weak coupling between the grains 
leads to a high thcrmopower and a low thermal conduc- 
tivity resulting in relatively high values of the figure of 
merit. We also discussed the case of compensated (half 
filling, v = 1/2) granular semiconductors. We showed 
that in this regime the thermopower is finite, but rather 
small. 

A. G. was supported by the U.S. Department of En- 
ergy Office of Science under the Contract No. DE-AC02- 
06CH11357. 
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